Effect of temperature on resonant electron transport through stochastic conduction 

channels in superlattices. 
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We show that resonant electron transport in semiconductor superlattices with an applied electric 
and tilted magnetic field can, surprisingly, become more pronounced as the lattice and conduction 
electron temperature increases from 4.2 K to room temperature and beyond. It has previously been 
demonstrated that at certain critical field parameters, the semiclassical trajectories of electrons in the 
lowest miniband of the superlattice change abruptly from fully localised to completely unbounded. 
The unbounded electron orbits propagate through intricate web patterns, known as stochastic webs, 
in phase space, which act as conduction channels for the electrons and produce a series of resonant 
peaks in the electron drift velocity versus electric field curves. Here, we show that increasing the 
lattice temperature strengthens these resonant peaks due to a subtle interplay between thermal 
population of the conduction channels and transport along them. This enhances both the electron 
drift velocity and the influence of the stochastic webs on the current-voltage characteristics, which we 
calculate by making self-consistent solutions of the coupled electron transport and Poisson equations 
throughout the superlattice. These solutions reveal that increasing the temperature also transforms 
the collective electron dynamics by changing both the threshold voltage required for the onset of 
self-sustained current oscillations, produced by propagating charge domains, and the oscillation 
frequency. 

PACS numbers: 73.21.-b, 05.45.Mt, 72.20.Ht 
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I. INTRODUCTION 



Semiconductor superlattices (SLs) are nanostructures 
formed from several alternating layers of different semi- 
conductor materials^"— . This periodic structure leads to 
the formation of energy minibands that enable electrons, 
in the presence of an electric field, to demonstrate a 
number of interesting quantum-mechanical phenomena, 
which include the formation of Wannier-Stark ladders, 
sequential and resonant tunneling, Bragg reflections, and 
Bloch oscillations. Consequently, SLs are of a great inter- 
est for both fundamental and applied scienceii^"— . Due 
to the high mobility of miniband electrons and the very 
high frequency of the Bloch and charge-domain oscilla- 
tions, SLs have prospective applications in sub-THz and 
THz electronic devices2ii£iiiii^iii. 

Recently, it has been shown that a tilted magnetic field 
applied to a SL can strongly affect, and hence control, 
the electrical properties of the structure. Nonlinear in- 
teraction between the electronic Bloch oscillations along 
the superlattice and cyclotron motion in the plane of the 
layers induces chaotic semiclassical electron dynamics, 
which, depending on the ratio between the Bloch and cy- 
clotron frequencies, either accelerate or decelerate charge 
transport through the S L^^d^ . Coupling of cyclotron and 
Bloch motion by a tilted magnetic field has also been 
related to the Fiske effect in superconducting Joseph- 
son junctions^S. On resonance, when the Bloch and cy- 
clotron frequencies are commensurate, the electrons ex- 



hibit a unique type of quantum chaos, which does not 
obey Kolmogorov-Arnold-Moser theory^^. This type of 
chaos is characterised by the formation of intricate "web- 
like" structures, known in the literature as "stochastic 
webs"—"—, which extend throughout the phase space 
of the miniband electrons. The appearance of these 
webs abruptly delocalises the electrons in real space, thus 
significantly increasing their drift velocity^^'^^. Simi- 
lar dynamics can occur in other spatially periodic sys- 
tems in which wave interference gives rise to band trans- 
port phenomena, including ultracold atoms in an opti- 
cal lattice^^, graphene^^, and light propagating through 
spatially- modulated photonic crystals^^'^^. On the quan- 
tum level, such resonant delocalisation of the electrons 
manifests itself in the formation of additional magnetic- 
field-induced miniband structure corresponding to ex- 
tended electron states, which extend across many periods 
of the SL^. This delocalisation dramatically affects the 
collective electron behavior by inducing multiple prop- 
agating charge domains and GHz-THz current oscilla- 
tions with frequencies and amplitudes much higher than 
with no tilted field^°. It has also been shown that in the 
vicinity of Bloch-cyclotron resonances the usually unsta- 
ble Bloch gain profile becomes stablc'^^, which can be 
used for the amplification of THz signals. Although the 
above effects of the Bloch-cyclotron resonances on elec- 
tron dynamics have been reported for different finite tem- 
peratures, it is still unclear how the character of charge 
transport in the SL changes with temperature. 

In this paper, we study how increasing the electron and 
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FIG. 1: (Color online) Schematic diagram of a semiconductor 
superlattice with an electric field and a tilted magnetic field 
applied. Red (green) arrows indicate co-ordinate axes (field 
orientation). Blue (yellow) bands in the y — z plane repre- 
sent quantum wells (tunnel barriers) within the superlattice. 
Electric field, F, is applied anti-parallel to the a;-axis and the 
magnetic field, B, lies in the (a;,2)-plane at an angle 6 to the 
a;-axis. 

lattice temperature affects tlie drift velocity of tlie elec- 
trons in an electric and tilted magnetic field and, conse- 
quently, tlie electric current through the SL. We find that 
the thermal distribution of the electrons can, counter- 
intuitively, enhance the effect of Bloch-cyclotron reso- 
nances on the drift velocity. In particular, increasing the 
temperature makes the resonant maxima in the drift ve- 
locity versus electric field characteristics more prominent. 
This shifts the frequency and amplitude of the collective 
charge-domain oscillations, associated with each resonant 
peak, together with the threshold voltage at which the 
oscillations switch on. 



II. ELECTRON DRIFT VELOCITY 

We consider the field configuration shown schemat- 
ically in Fig. [H The electric field F = (-F,0,0) 
is applied perpendicular to the plane of the SL lay- 
ers and anti-parallel to the x-axis. The magnetic field 
B — {B cos 9,0, B sin 9) lies in the (x, z)-plane at an an- 
gle 9 to the a;-axis. The semiclassical equation of motion 
for a miniband electron is 

p(<) = -e[F + (Vpi?(p(t)) xB)] (1) 

where e is the magnitude of the electron charge, p{t) — 
{px{t) , Py{t) , pz{t)) is the electron's crystal momentum at 
time t and E{p) = A{1 - cos{p^d/hj)/2 + (p^ + pi) /2m* 
is the dispersion relation for the lowest miniband within 
the tight-binding approximation. Here, A is the mini- 
band width, d is the SL period, and m* is the electronic 
effective mass for motion in the (y, z)-plane. 



Equation ([T]) can be expressed in its constituent 
componentsi^i^ 

p,(t) = eF-uj^py{t) (2) 

. dAm*ujj_ . (p^{t)d\ 
Pyi*) = 2l sm ( — ^ I - W||p^(i) (3) 

Mt) - ^WPyit), (4) 

where LO\\=eB cos 9/111,* and uj±=eB sin 9/ m* are the cy- 
clotron frequencies corresponding to the magnetic field 
components along the x- and z-axes, respectively. The 
electron velocity along the x-direction is given by 

v^(t)=i{t)=vosin(^P^y (5) 

where the peak velocity vq = Ad/{2h). 

To determine the drift velocity, Ud, of electrons with 
initial momentum P — {Px, Py,Pz) we use the Esaki-Tsu 
formalisntiii^: 

POO 

Ud{P) = iy vx{t)e-'''dt, (6) 
Jo 

where 1/ is the electron scattering rate^. 

In the case of non-zero temperature, T, one should 
take into account the thermal distribution of electron 
momenta f (P^ , Py , Pz) , which we assume to obey the 
Boltzmann statistics^i^i^i 

f(P'j^±(,~2k^[^-''°^—)-27^Pn^_ (7) 

In Eq. ([7]), Z can be found from the normalisation con- 
dition 

/Trh/d />oo poc 
/ / f{V)dPxdPydP, = 1, 

-Trh/d J —00 J — 00 

which yields 

Z = {2T:fm*kBT-^h (^) e"^, (8) 

where Iq{x) is a modified Bessel function of the first kind. 
Then, the drift velocity of miniband electrons at temper- 
ature T is 

/ Trh/d />oo 1*00 
/ / f{P)Ud{P)dPxdPydP,. (9) 

-Trh/d J —00 J —00 

In the Appendix we show, following previous work^^^, 
that for small angles 9, for which a;j_ ^ wy and uj^ <g; 
cjs, where lob — eFd/h is the frequency of the Bloch 
oscillations, the drift velocity can be approximated by 
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Ji(A/2fcBT) 
/o(A/2fcsT) 



TO kBi — -T 

Vw|| ft, 



v{ujb — nLj\\ ) 
' + [lob - nuj\\y 



(10) 



where /„ (x), = 1, 2 . . . are the modified Bessel functions 
of the first kind. When B — Q, Eq. ([TU|) reduces to 



/i(A/2fcBr) t/cjB 
Io{A/2kBT) iy^+uj\ 



(11) 



Note that Eq. ([TT]) can be derived exactly^ by inte- 
grating Eqs. ©-(HI) with 5 = 0, and substituting the 
resulting expressions for p^, Py and p^ into Eqs. ([5]) and 
^ . For T = the equation ([TT|) yields the famous Esaki- 
Tsu relation for the electron drift velocity^. 

In the general case of arbitrary parameters F, B and 
6', Eqs. ©-(HI, which, in principle, can exhibit deter- 
ministic chaoai^ii^i^, cannot be solved analytically, and 
therefore the electron drift velocity must be calculated 
numerically. In our computations we use the following SL 
parameters, taken from recent experiment a^^'^^ : d = 8.3 
nm, A = 19.1 meV, = 4 x 10^^ s^^ and n 



0.067top 



where 



III. 



is the mass of a free electron. 



NUMERICAL CALCULATION OF THE 
ELECTRON DRIFT VELOCITY 



In our numerical simulations, we integrate Eqs. ([2])- 
(|3]) for a number of different initial conditions p(0) = P, 
which were randomly generated according to the proba- 
bility density function ([7]). This function can be decom- 
posed into momentum components 



/(P) = 



1 



^/27^TO*fc^ 
1 

y/2Trm*kBT 
d 



exp 



exp 



2m*kBT 
2m*kBT ^ 



(12) 



2TThIo 
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\2kT 



A fP^d 



As a consequence of this decomposition, one can use 
independent random number generators for each of the 
initial momentum components P^, Py, Pz- The values 
of Py and Pz can be obtained from standard routines^ 
for the generation of uncorrelated random numbers obey- 
ing a Gaussian distribution with zero mean and variance 
(fc^TTO*)^/^. However, the generation of random values 
of Px requires a more sophisticated procedure, which we 
now explain. 

Let us find a continuous function, g{S,), which trans- 
forms a random variable, ^, uniformly distributed within 



the interval 
distribution 



-7r,7r), into a random variable, g, having a 



dexp 



/(5) = 



2kBT 



27rhIo I I 
\2kBTj 



exp 



A 



2kBT 



1 I Sd 

1 — cos I — 
h 



(13) 

To find 5(0, we use the probability conservation re- 
lation (l/27r)(i^ — f{g)dg, which gives us a differential 
equation for g{^): 
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(14) 



dgiO 



Due to the symmetry of the distribution in Eq. (jTSj), 
g{0) = 0, which can be used as an initial condition for Eq. 

Moreover, it is possible to show that (7(±7r) — zLtt 
for any values of temperature T. However, in the general 
case of arbitrary argument, f , g{^) cannot be expressed 
analytically. We therefore evaluate it by numerical in- 
tegration of Eq. Figure m illustrates the function 
g{^) for three different temperatures T = 50 K (curve 1), 
T = 150 K (curve 2), T = 300 K (curve 3). As the figure 
shows, g{^) is an odd function, which becomes almost 
linear as T increases. 

In practice, numerical integration of Eq. (1141) signifi- 
cantly slows the simulations. Therefore, we used approx- 
imate analytical solutions. We found that for small tem- 
peratures, T < 50 K, the distribution can be well 
approximated by a Gaussian with zero mean and vari- 
ance 2kBTh/{Ad). Hence, in this case P^, can be also 
defined using the Gaussian random number generator. 
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FIG. 2: (Color online) Function g{£^) for three different tem- 
peratures: T = 50 K (curve 1), T = 150 K (curve 2), T = 300 
K (curve 3). d = 8.3 nm, A = 19.1 meV, : 
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TABLE I: Coefficients for tlie polynomial expansion in Eq. 
(|15|l at different temperatures T 



r 


ai 


13 


as 


a? 


WOK 


0.42008 


0.04124 


-0.00755 


0.00094 


125K 


0.49513 


0.02223 


0.00087 


0.00039 


150K 


0.55066 


0.01556 


0.00210 


0.00010 


175K 


0.59097 


0.01440 


0.00309 


-0.00003 


200K 


0.62431 


0.01505 


0.00329 


-0.00009 


225K 


0.65325 


0.01626 


0.00312 


-0.00012 


250K 


0.67807 


0.01744 


0.00282 


-0.00013 


275K 


0.69868 


0.01836 


0.00249 


-0.00013 


300A' 


0.71677 


0.01905 


0.00217 


-0.00012 


325K 


0.73367 


0.01954 


0.00186 


-0.00011 


350K 


0.74852 


0.01983 


0.00159 


-0.00010 


375K 


0.76150 


0.01996 


0.00136 


-0.00009 


400 A' 


0.77299 


0.01998 


0.00116 


-0.00009 



For larger temperatures T > 50 K we approximate the 
solution of Eq. by the following polynomial expres- 
sion 



/(O = aiC + + + ^7^' + Oif 



(15) 



where the coefficients a„ for different temperatures are 
given m Tab. n Note that with increasing T, the coeffi- 
cient fli, corresponding to the linear term, grows, whereas 
the other (nonlinear) coefficients, 03, 05, 07, decrease. 
Also, the higher-order terms can be neglected. 

Thus, to calculate the drift velocity, Vd, for given tem- 
perature, T, we (i) randomly generate the components 
of the initial momentum vector P, (ii) use them as an 
initial condition to integrate Eqs. (O-Q, (iii) substitute 
the results of this integration into Eqs. ([5]) and ([6]) to 
determine the drift velocity, Vd, for the given initial mo- 
mentum, and (iv) average Vd over all randomly-chosen 
initial momenta. 

To ensure that the above procedure for simulating 
the thermal distribution of electron momenta is accurate 
enough, we check it for the case of no applied magnetic 
field, for which Vd can be calculated analytically^ using 
Eq. pT|) . The solid curves in Fig. ^a) shows Vd{F) plots 
calculated using Eq. pT|) for T = 4.2 K (upper curve) 
and T = 300 K (lower curve). The data points marked by 
symbols "•" and "o" show values of Vd obtained numeri- 
cally for T = 4.2 K and T = 300 K, respectively. The fig- 
ure demonstrates excellent agreement between the exact 
Eq. (jlll) and our numerical procedure. More precisely, 
we find that by averaging over N > 200000 trajectories, 
the relative error does not exceed 1 — 2%. Consequently, 
in all our simulations we used N = 250000. 

Eq. pT|) shows that the maximum drift velocity, ti™"^ ^ 
is achieved when tug = v, i.e. when a significant fraction 
of the electrons complete whole Bloch oscillations before 
scattering. As lub increases beyond v, v^""^ decreases 
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FIG. 3: (Color online) (a) Variation of drift velocity, Vd, with 
electric field, F, determined analytically from Eq. (Ilip (solid 
curves) and calculated numericaUy (symbols) for T — 4.2 K 
(curve 1 and "•" ) and T = 300 K (curve 2 and "o"). (b) 
Variation of maximal drift velocity, w™"^, with T calculated 
analytically (sofid line) and numerically (•). Numerical sim- 
ulations involved averaging over 250000 different initial mo- 
menta, P, chosen randomly. 



with increasing T according to the equation 

A 



2kBT 



A 



2kBT 



(16) 



The variation of with T predicted by Eq. (|T6)) 

is shown by the solid curve in Fig. |31Ib). For compari- 
son, the result of our numerical simulation is shown by 
the symbols "•" . The figure shows excellent quantitative 
agreement between our analytical and numerical results. 



IV. EFFECT OF TEMPERATURE ON 
ELECTRON DRIFT VELOCITY IN A TILTED 
MAGNETIC FIELD 

We have also applied the numerical procedure de- 
scribed in the previous section to the case of a tilted mag- 
netic field, for which the drift velocity Vd cannot be de- 
termined analytically. Figure 2] shows Vd versus F curves 
calculated numerically for B = 15 T and 9 — 40° at 
several different temperatures T. 

In contrast to the case of i? = [Fig. EJa)], all of 
the Vd{F) curves in Fig. [4]exhibit multiple maxima. The 
first maximum, for the lowest value of F, also exists when 
B = and, as noted in the previous section, is associated 
with the onset of Bloch oscillations. All other maxima oc- 
cur because of the enhanced acceleration of the electrons 
whenever the ratio of the Bloch and cyclotron frequencies 
r = ujb/^± = 0.5, 1, 2, or 3 (upper scale in Fig. 
Remarkably, as T increases, the first peak weakens signif- 
icantly, whereas the amplitude of the cyclotron-Bloch res- 
onances (i.e. the peak-to-valley ratio) increases sharply 
even though the peak Vd values decrease. Moreover, for 
high enough temperatures, new peaks, reflecting higher 
order resonances, appear in the Vd{F) dependencies. 

To qualitatively understand the effects of tempera- 
ture on drift velocity, we analysed single electron or- 
bits described by the Eqs. ©-([l]). Figure [51^ a) shows 
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FIG. 4: (Color online) Variation of drift velocity, Vd, with 
(lower scale) electric field, F, (upper scale) r = ujb/oj\\ calcu- 
lated for B = 15 T and 9 = 40° and different temperatures. 
Curve 1 corresponds to T = K; curve 2 to T = 50 K; curve 3 
to T = 200 K; curve 4 to T = 300 K; curve 5 to T = 400 K. 



a stroboscopic Poincare section phase portrait for the 
r = LOB/ijJ\\ = 3 resonance (S = 15 T and 9 ~ 40°), 
taking the strobe period to be Tb = 2Tr /lob- In ad- 
dition to localised quasi-periodic orbits, which appear 
in the stroboscopic section as dotted closed curves, the 
phase space also contains a resonant unbounded structure 
formed by chaotic layers and known in the literature as a 
"stochastic web"— It has been shown previousl y 
that electrons having initial momenta in the stochastic 
web move rapidly through the SL, so producing a large 
resonant peak in the Vd{F) curve, as shown in Fig. |4] 
for r = 0.5,1,2,3,4. The height of each resonant peak 
increases as more electrons have momenta lying within 
the stochastic web. In Figs. [SJb),(c) we show enlarged 
fragments of the stochastic web near the origin together 
with a grey-scale map indicating the probability of find- 
ing an electron with given initial momentum components 
iPy,Pz) at T = 50 K [Fig. [5] (b)] and T = 400 K [Fig. 
[SJc)]. Comparison of these two figures reveals that, with 
increasing temperature, the probability of the electron's 
initial momentum lying within the stochastic web grows. 
At the same time, the number of electrons with higher 
initial momenta, and thus higher velocity, also increases. 
Together, these two effects increase Vd and thus enhance 
the Bloch - cyclotron resonant peaks in Fig. 

Interestingly, although Eq. (ITOl) is unable correctly to 
predict the exact electron drift velocity for an arbitrary 
magnetic field configuration, it still captures certain qual- 
itative trends in the variation of Vd with the field param- 
eters and T. First, we note that according to Eq. (fTOt . 
the expression for Vd consists of a series of resonant terms 
proportional toiy{u!B ±«W| | ) / [h''^+{ujB^nu!\ | )^] . The term 
with n — reflects the Esaki-Tsu Vd{F) curve, whereas 
other terms with n ^ correspond to the Bloch - cy- 
clotron resonances, for which ujb/^^w = n is an integer. 
Moreover, since, for very large x, In{x) has the asymp- 
totic form In{x) oc exp(a;)/v^, as T increases in Eq. PH)) 
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FIG. 5: Stroboscopic Poincare section of electron trajectories 
in the {Py,Pz) plane taken with strobe period Tb ~ 2tt/lub 
for the r = a;_B/aj|| — 3 resonance {B — 15 T, 6 = 40°). (a) 
"Stochastic web" and neighboring quasi-periodic orbits; (b), 
(c) enlarged part of "stochastic web" (black dots) and prob- 
ability (grey-scale map in which lighter shades correspond to 
higher probability) that the initial momentum {Py, P^) lies in 
the given area of the {py,Pz) plane calculated for T = 50 K 
(b), and T = 400 K (c). 



the higher-order resonant terms strengthen. Under cer- 
tain conditions this leads to the appearance of additional 
peaks in Vd{F), which we see in our numerical simulations 
shown in Fig. [H 
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V. EFFECT OF TEMPERATURE ON 
ELECTRIC CURRENT THROUGH THE 
SEMICONDUCTOR SUPERLATTICE. 

In order to study the effects of temperature on the 
collective dynamics of electrons in the SL, we self- 
consistently solve discrete versions of the current con- 
tinuity and Poisson equations, splitting the miniband 
transport region into N — 480 layers, each of width 
Aa; = L/N = 0.24 nm, which is small enough to ap- 
proximate a continuun>2fl"i^ . 

The evolution of the electron density, Um, in the m*'' 
layer, whose right-hand edge is at a; = mAx, is given by 
the discretised current continuity equation 

eAx—^^Jm-i-Jm, m = l...iV, (17) 

where Jm-i [Jm] are the areal current densities at the 
left- [right-] hand edges of the m*'' layer. Within the 
drift-diffusion approximation, the current density is 



erimVd 



Ax 



(18) 



where Fm is the mean electric field in the m*^ layer- 
ed {Fm) d ( eFmd 



D (Fm) = 



1 — exp 



( eFmd \ 



exp 



kT 



(19) 



is the diffusion coefficient'^, and the drift velocity {Frn) 
is obtained by using the approach described in Section 
mil Since Jm depends on the local drift velocity, Vd (Fm) , 
the collective electron dynamics depend directly on the 
single-electron orbits. 

The electric field Fm [Fm+i] at the left- [right-] hand 
edges of the to*'' layer can be described by the discretised 
Poisson equation 



F 



m+l 



eAx 

eo£r 



{um - no) + F„ 



I...N, (20) 



where Eq and £r — 12.5 are, respectively, the absolute 
and relative permittivities and no = 3 x 10^^ m~^ is the 
n-type doping density in the SL layers^^. In the emitter 
and collector Ohmic contacts, F — Fq. 

We use Ohmic boundary conditions to determine the 
current, Jo = ctFq, in the heavily doped emitter contact 
whose electrical conductivity a = 3788 Sm~^.^^ The volt- 
age, V, applied to the device is a global constraint given 
by 



Ax 
'2 



N 



V=U+— y(Fm+Fm+l), 



(21) 



where the voltage, U, dropped across the contacts in- 
cludes the effect of charge accumulation and depletion in 
the emitter and collector regions and a contact resistance 
i? = 17 Jl^. We calculate the current as 



I{t) 



N+1 



N 

E 



(22) 
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FIG. 6: (Color online) I{V) characteristics calculated for 
B = T (a) and B = 15 T, e = 40° (b) at (from top to bot- 
tom) T = K, 50 K, 100 K, 200 K and 300 K. In (a), data 
points on dashed line show values of {Vth, hh) calculated from 
Eqs. (|23|l and (|26p and discussed in the text. 



where A = 5x10 m^ is the cross-sectional area of the 
SL.^^^ 

The model described by Eqs. (HIl) - (l22|) exhibits both 
constant and oscillating electric current, depending on 
the voltage, V, applied to the device. Typical current- 
voltage [/(l^)] characteristics are shown in Fig. |6] both 
for i? and when a tilted magnetic field {B — 15 T, 
6 = 40°) is applied. Note that for V values at which 
current oscillations occur, the DC current was calculated 
by averaging I{t) over time. 

When B = [Fig. [5] (a)], the I{V) curves reveal the 
usual Esaki-Tsu-like behavior, characterised by a single 
maximum, which is associated with the onset of single- 
electron Bloch oscillations throughout much of the SL 
transport region. Figure [6] (a) shows that as T increases, 
the peak current decreases, which agrees well with a num- 
ber of experimental observations'^^''**'. This decrease of 
the peak current reflects that of the maximal electron 
drift velocity, predicted by Eq. ([T6l) . as T increases. In 
particular, the factor h{A/2kBT)/ Io{A/2kBT), charac- 
terising the effect of the temperature on the drift veloc- 
ity in Eq. (fT6|). also accurately describes the drop of 
the peak current both in our simulations and in earlier 
experiments discussed in Refs<^i^ 



The electron dynamics changes significantly when a 
tilted magnetic field is applied to the SL. Figure [5] (b) 
shows typical I{V) cmves calculated, in the presence of 
a tilted magnetic field, for a range of temperatures. All of 
the I{V) curves reveal clear Bloch-cyclotron resonances, 
which manifest themselves through the appearance of ad- 
ditional features in the curves^. For low temperatures, 
the Bloch-cyclotron resonances produce sudden changes 
in the slope of the I{V) curves [for example, the kink 
near V = 0.3 V in the (upper) I{V) curve calculated for 
T = in Figure [6] (b)] and also shift the position of the 
current peak. As T increases, these effects become more 
prominent and, eventually, give rise to additional maxima 
in the I{V) curves, e.g. the double peaks at T = 200 K. 
This evolution of the I{V) characteristics originates from 
the variation of the Vd{F) curves with changing temper- 
ature, shown in Fig. |4l Comparison of Figs. [6] (a) and 
(b) shows that although the maximum current decreases 
with increasing temperature in both cases, the peak DC 
current is always larger in the presence of a tilted mag- 
netic field. This can be explained by the strong resonant 
enhancement of Vd produced by the Bloch-cyclotron res- 
onances. Moreover, as Fig. 2] indicates, with increasing 
temperature the amplitude of these resonances increases 
and so the difference between the peak currents for i? = 
and for i? 7^ becomes larger. 

When the applied voltage, V, exceeds some critical 
value, Vth, which depends on T, B, and 9, the stationary 
state of the system loses its stability, and the electric 
current starts to oscillate at a frequency, /, which, for 
the given parameters, is in the microwave range. These 
current oscillations are associated with the formation of 
traveling charge domains'^''^^. 

The variation of the oscillation frequency, /, with V 
is shown in Fig. [7] for a range of T. Note that for each 
T value, the corresponding f{V) curve starts at a par- 
ticular threshold voltage, Vth, above which /(t) oscilla- 
tions occur. In general, as V increases, / decreases for 
both 5 = [Fig. El (a)] and B ^ 15 T, 9 = 40° [Fig. [7] 
(b)]. This can be understood in terms of the Vd{F) curves 
shown in Figs. OandlH These figures suggest that when 
F is sufficiently large, Vd decreases with increasing F for 
both B — and i? = 15 T. Moreover, further increas- 
ing F in this regime slows the electrons even more. The 
applied voltage produces an inhomogeneous electric field 
distribution throughout the SL layers, which is deter- 
mined by Eqs. p?l) - ((2T|) . As V increases, the region of 
high electric field inside the SL, where Vd is small, also 
increases and so the speed of the charge domains is low. 
The larger the region of high electric field, the smaller the 
average velocity of the charge domains and, hence, the 
lower the frequency of current oscillations. When B = 15 
T and 9 = 40°, the multiple peaks in the Vd{F) curve 
tend to keep the average Vd value high, which means 
that the charge domains move faster than for B ^ 0. 
In addition, the onset of Bloch-cyclotron resonances in- 
creases the maximal Vd value that the electrons attain 
for given T [see Fig. 2] . All of these factors improve the 
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5=15 T 




F(V) 

FIG. 7: (Color online) Variation of the generation frequency, 
/, with voltage, V, applied to the SL calculated for _B = T 
(a) and B = 15 T, 6 = 40° (b) at the temperatures indicated. 



propagation speed of charge domains in the presence of 
a tilted magnetic field. Thus, for given V, the frequency 
of charge domain oscillations (if they occur at that V) is 
higher when a tilted magnetic field is applied. 

The evolution of the amplitude of the current oscilla- 
tions, A/, with V is shown in Fig. [51 Comparison of Figs. 
H] (a) and[8](b), reveals that applying a tilted magnetic 
field at any temperature allows one to generate current 
oscillations of higher amplitude and power. The mech- 
anisms that produce this power enhancement were dis- 
cussed in detail in Refi^. In particular, it was shown 
that multiple peaks in the Vd{F) curves induce multiple 
propagating charge domains, which strengthen the cur- 
rent oscillations. Since the multi-peak character of the 
Vd{F) curves persists over a wide range of T (see Fig. 
S]) , the power enhancement induced by a tilted magnetic 
field occurs even at room temperature. 

Remarkably, the modification of the I{V) curves by a 
tilted magnetic field can lead to the appearance of new 
regions of negative differential conductance, for example 
associated with the two peaks in the I{V) curve at T = 
200 K (purple curve in Fig. (5]). In this case, there are 
several distinct regions of V within which the SL exhibits 
current oscillations [see Fig. [S](b)]. 
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FIG. 8: (Color online) Variation of the amplitude of the cur- 
rent oscillation, A/, with voltage, V, applied to the SL cal- 
culated for B = T (a) and B ^ 15 T, 9 = 40° (b) at the 
temperatures indicated. 



Figures [7] and [5] indicate that for both B = and for 
_B = 15 T, the threshold voltage for current oscillations, 
Vth , depends non-trivially on T. For the low temperature 
range T < 100 K, increasing T decreases Vth- But as T 
increases further, Vth starts to increase. This is due to 
the effect of temperature on the current density, J, which 
is described by Eqs. ([TS]) and (fT9|). As noted above, the 
generation of oscillating current in the SL is associated 
with the formation of traveling change domains, which 
occurs when the voltage, Uf/^, dropped across the active 
layers of the SL becomes sufficiently larg o^'^° . Eq. (PT|) 
shows that this voltage is the difference between the total 
voltage, V, applied to the SL and the voltage, U, dropped 
across the contacts. The voltage, U, in turn, may be ex- 
pressed as [/ = IR + Uc, where Uc describes the effect of 
charge accumulation and depletion in the emitter and col- 
lector regions. Thus, any changes in the maximal value 
of J and, therefore, of /, will affect U, and thus shift the 
threshold voltage Vth- For low temperatures, when elec- 
tron diffusion can be neglected, increasing T decreases 
the maximal electron drift velocity and, thus, the max- 
imal value of I. This reduces the voltage, U, dropped 
across the contacts, since this is proportional to the cur- 
rent through those contacts. When U decreases but V 
remains constant, the voltage, Uf^, dropped across the 



active layers 
perature ran; 

The value 
current Ith = 
T and when 
and uses the 
current, Ith 



of the SL increases. Thus, within some tem- 
ge, increasing T decreases Vth- 
of Vth and the corresponding critical DC 
liVth) can be estimated analytically at low 
B = 0. If one neglects electron diffusion 
drift velocity given by Eq. (fT6|) . the critical 
at temperature T can be estimated as 



Ith{T)=I?^ 



h{A/2kBT) 
/o(A/2fcsT)' 



(23) 



where I^j^ is the critical current for T = K. According to 
Eq. ([2T|) . the critical voltage, VthiT), required for current 
oscillations is 



Vth{T)^U, + IthiT)R + U, 



SL 
th J 



where 



SL 
th 



Ax 
~2~ 



N 



(24) 



(25) 



m—l 



is the voltage dropped across the SL layers. Assuming 

tSL 
^th 



that Uc and UfJ" are almost independent of T, one ob- 



tains 



VthiT) V,, - I,,R \l ^^(^/2fcsT) ) ' 



(26) 



where is the critical voltage at T = 0. From Eqs. 
(I23t and dMl) it follows that 



i?(/,.(T)-/°J, {Ith{T)<ll). (27) 



Vth{T) 



V^ 



th 



The Ith{Vth) curve defined by Eq. (f27| is shown by 
the dashed line in Fig. |6] (a). The open circles mark 
the coordinates {Vth{T) Jth{T)) for the five different T 
values corresponding to the five I{V) plots shown in the 
figure. The values of and were taken from the 
I{V) curve for T = 0. From Figs. H (a) andlTia) it is 
easy to see that for the low T values considered (< 300 
K), the analytical values of Vth and Ith obtained from Eq. 
(|27p are in excellent quantitative agreement with those 
inferred from the numerical I{V) curves. 

For large T, the diffusion contribution to the current 
becomes important and Eq. p7| fails. As T increases, 
the current through the SL and the voltage dropped 
across the contacts also increase, thereby raising Vth- 

The same mechanism operates when _B 7^ 0. In this 
case, Vth is larger than for _B = (see Fig. [T]). This 
is because in the presence of a tilted magnetic field, the 
maximal drift velocity is larger due to the onset of Bloch- 
cyclotron resonances (compare Figs. 3 and 4). 



VI. CONCLUSIONS 

In conclusion, we have shown that a tilted magnetic 
field strongly affects, and can significantly enhance, the 
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transport characteristics of SLs - even at room temper- 
ature. The temperature dependence of the electron drift 
velocity in the presence of a tilted magnetic field is quite 
different from the B — case. In particular, increasing 
T quickly suppresses the B — Esaki-Tsu peak in the 
Vd{F) curve but has a much smaller effect when a tilted 
magnetic field is applied. In this geometry, increasing T 
sometimes even enhances the drift velocity peaks, relative 
to the background, caused by the Bloch-cyclotron reso- 
nances. As a result, when B and 9 are both non-zero, the 
electrons move faster, which increases both the DC- and 
AC-components of the current through the SL and also 
improves its high-frequency performance. The effects of 
temperature on the Vd{F) curves also transform the SL's 
I{V) curves by inducing new instabilities and transport 
regimes, which we hope will stimulate future experimen- 
tal and theoretical studies. Since the effects that we have 
demonstrated are generic features of semiclassical energy 
band transport, we expect similar phenomena to occur 
in, for example, cold atoms^^'^^ and photonic crystals^!. 
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VIII. APPENDIX 

We solved Eqs. ©-(Hj) assuming a;|| Under this 

condition Eqs. ©-(HI) can be linearised as 

p,{t) = eF-Lj^py (28) 
Pyit) ^ ~LUi\p,{t) (29) 
p,{t) = uj^^pyit), (30) 

which for initial conditions Pi(0) = p||, Py{0) = Py, 
p,{0) = yields 

Px{t) =P\\ +eFt p±_ sm(u}Ht + tpo) -\ pi_ sm(^o, 

Pyit) ^ p±cos{u!\\t + ipo), (31) 
Pzit) ^ p±sm{uj\\t + ipo), 



where p± = yPy + Pzi a-nd Lp^ — atan{Pz/ Py). 

Eq. (j31l) for the crystal momentum can be used to 
calculate the electron velocity along the SL (x) axis as a 
function of time, t: 



Vx [t) = Wo sin 



Px{t)d 



Wo 
2i 



P\\d uji_ d LJ±_ d 

vo sm — — + LUBt -pi_ sm(a;||t + ipa) H -pi_ sm(y9o 

n n cj|| h 

I P\\d UJ± d Lu± d 

exp I—— + iLOst — « — tPI- sm(cji|t + csq) + « — TPi- smcso 
' n. uj\\ n " W|| n 

I P\\d Ll>± d Lu± d 

exp —I— lujBt + i — tP± sm(uj»t + ipo) — i — -p± smipo 

' n. a;|| n. " W|| n 



(32) 



r 



where ujb = eFd/h is frequency of the Bloch oscillations Using the Jacobi- Anger expansion for the Bcsscl func- 

and vq — Ad/{2h) is the maximal velocity in the given tions we obtain 

miniband. 



Mt) = ^^ 



n— — ooA:— — oo ^ ii ^ ii 



i{pQ{n — k) 



(33) 



n= — cxD/c= — oo ^ I 



r 



The drift velocity, Ud(jp\\,p±), of electrons with initial mo- mentum components p\\ and p_L can be calculated using 
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the Esaki-Tsu formula (|6]). Substituting (l33l) in ^ we get 



2,1 



P||d 



(34) 



with 



oo oo 



n— — oo A;— — c 



r 



To find the drift velocity, Vd, at finite temperature T we tion of electron momenta ([7]). In the current notations 
average Ud{p\\ ,P±,, V?o) according to the thermal distribu- the definition of Vd becomes 



Vd 



oo rTT 



ttH/ d 



irh/ d J —OQ J —t: 



Ud{p\\,p±^,(po)fiP\\,PJL)p±dpi\dpj_dipo, 



(35) 



r 



with 



f{p\\,Pi-) = ^exp 



A / P\\d 

1 — cos 

2kBT \ h 



Pi 



2m*kBT 



To determine Vd from the above integral expression, we 
first evaluate the integral over p|[, which can be written 
as 



and Z is given by Eq.(| 



or 



1 vo A_ r^/'^ 

J-TTh/d 



UdiP±,Vo) 



Z 2i 
' Z2i 

1 Vo 



p\\d' 
exp [ I'-j— ] exp 



Tzh/ d 

B{p±,ipQ)e~^^ I exp 

7rh/d 



pTTh/d / A \ 

J-^n/d^t^^ \2kBTj 

/Trh/d oo / A \ 



A cos 



Piid 



.P\\d' 
-^— ) exp 



2kBT 
"Acos^ 



2kBT 



dP\\- 



77 /l/ d 



Z 2i 



2-K Vq h _ 



/ A 



\2kBT 
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At the next step we perform integration over p± and 



Vd 



Ud{p±,(po) exp 



J -TT 



Pa 



2m*kBT 



p±dp±d(po 



27r vq h A / A 



n— — oo fc^ — oo 





r ? n 

PI 


exp 


2m*kBT 



(27r)2wo n 



-e 



p±dp±dipa = 



X exp 



9 

PI 



2m*kBT 



p±dpA 



r 



(36) 



Taking into account that J-n{x) ~ (— l)"J„(a;), we can rewrite (|36 
I 



n=0 



^" ' —^P^ ) -,2^f,^,^_„ ^2 X exp 



uj±_ d 

— P 

uin n 



2 fujj_ d 

Jn TP^ 



v'^ + {ujB — nijj\\)- 



P\ 



I 



2 X exp 



2m*kBT 

„2 



2m*kBT 



p±dp±_ + (37) 



The integral 



p_L exp 



Pa 



2m*kBT 



P±dpj 



I 



-m*kBT 



uj± d 

cjii h 



m*kBT 



uj±_ d 

LOW ^ 



(38) 



is, for n > —1, a special from of the second Weber ex- 
ponential integral^. After integration and re-indexing 



J 



the sum in (jST]), and substitution of the expression for Z 
from we obtain the final formula for the drift velocity 



/i(A/2fcsT) 
^^ = "Vo(A/2fc.T)^^P 



-m kBi 



n— — 00 



m kBi 

W|| h 



v'^ + {lob - JT-i^ii)^ 
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